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D . An 50(10) SUSY GUT model was previously presented based on a minimal set of 
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Higgs fields. The quark and lepton mass matrices derived fitted the data extremely well 
and led to large f^ — i^r mixing in agreement with the atmospheric neutrino data and 
to the small- angle MSW solution for the solar neutrinos. Here we show how a slight 
modification leading to a non-zero up quark mass can result in bimaximal mixing for 
the atmospheric and solar neutrinos. The "just-so" vacuum solution is slightly favored 
over the large-angle MSW solution on the basis of the hierarchy required for the right- 
handed Majorana matrix and the more nearly-maximal mixing angles obtained. 
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Atmospheric neutrino data ||I|] reveal a large and, in fact, nearly maximal mixing between 
z/^ and some other neutrino, which most plausibly is assumed to be Ut-, though it may also 
be a sterile neutrino. The solar neutrino problem |0, on the other hand, can be solved either 
by a small mixing (sin^2^e^ ~ 6 x 10^^) or by a nearly maximal mixing of z/g with z/^. In 
the former case one has the small-angle MSW solution 0, while in the latter case one has 
either the large-angle MSW solution [^] or the vacuum oscillation solution [^. Cases where 
both the atmospheric and solar neutrino anomalies are solved by nearly maximal mixing |Q 
are often called "bimaximal." 

There have been several interesting suggestions about how bimaximal mixing might arise 
theoretically 0. Here we suggest a scheme that has certain novel features. The basic idea, 
which at first sounds artificial but which actually can emerge quite simply and naturally 
as will be seen, is that the large mixing of u^ and u^ originates from transformation of the 
charged lepton mass matrix, whereas the large mixing of z/g and z/^ originates in the neutrino 
mass matrix itself. (Of course, it only makes sense to draw this distinction if there is a 
preferred basis of families, which is here provided by some underlying theory of flavor.) The 
idea that the large z/^ — Ur mixing arises from the mass matrix of the charged leptons was 
proposed in 0, where it emerged as part of a complete model for the quark and lepton 
masses and mixings. One of the virtues of this idea is that it allows a simple resolution 
of the supposed paradox that the mixing of the second and third families is small for the 
quarks, Vcb — 0.04, and large for the leptons, V^s ~ 0.7. 

The model developed in [^ was only a model for the heaviest two families with the first 
family being approximated as massless. In [Q this model was extended to the first family 
in a very economical way that gave several additional predictions, among which were that 
the Uf. — Ufj_ mixing angle is small, and in fact precisely in the presently allowed range for the 
small-angle MSW solution. This small Ve ~ i^p. angle, like the large z/^ — Ur angle, arose from 
diagonalization of the charged lepton mass matrix. 

In the present paper we show that a slight alteration of that model leads to an equally 
predictive scheme that has bimaximal mixing. All the predictions for quark and charged 
lepton masses, for the CKM parameters, and for the z/^ — u^ mixing are left essentially 
unaffected; however, the Ug — v^ mixing moves from the small-angle MSW value to become 
maximal. We will first very briefly review the model of ^ and p], noting the features 
relevant to lepton mixing, and then proceed to show how bimaximal mixing can arise in it. 

The model is based on supersymmetric 50 (10) and leads to the following matrices: U^ 
for up-type quarks, D*^ for down-type quarks, V^ for the charged leptons, and N^ for the 
Dirac neutrino masses, where the superscript refers to the matrices at the unification scale. 
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These matrices arise from simple diagrams in the SO{10) unified model, which involve only 
five effective Yukawa terms that have the forms (163163)10/^-, (162163)10^1/45/^, [I62I6//] 
[16316'//], [16il62][16//16^], and [16il63][16//16'//]. These lead, respectively, to the entries 
in Eq. (1) that are denoted 1, 6, a, 6, and 6'. The numerical subscripts are family indices, 
while the subscript H denotes a Higgs multiplet. The notation [16 16] implies the spinors 
are contracted into an 5*0(10) vector. The VEV, (16//) ~ Mq, lies in the SU{5) singlet 
direction and helps break SO{10) down to the Standard Model, while (16^) ~ Mw lies in 
the weak doublet direction and helps break the electroweak interactions. The expectation 
value of the adjoint ((45//)) is proportional to the 50(10) generator B—L, as required for the 
Dimopoulos-Wilczek mechanism to provide the doublet-triplet splitting |p. The foregoing 
information is sufficient to derive the matrices in Eq. (1). 

The family hierarchy results from the smallness of the parameters 6 = 0.14, 6 = 0.008, 
and \6'\ = 0.008. The parameter a = 1.8 is not small, however, and is the key to under- 
standing many of the qualitative and quantitative features of the quark and lepton spectrum, 
including the fact that the second and third families of leptons have a large mixing while 
that of the quarks is small. The point is that SU{5) relates left-handed (right-handed) down 
quarks to right-handed (left-handed) charged leptons, and consequently relates the ij ele- 
ment of D^ to the ji element of L°. That is why in Eq. (1) the large parameter a appears 
in L32 where it leads to large /i~ — r~ mixing and hence large sin^26'^^, whereas it does 
not appear in D32 where it would give large Vet, but rather appears in D23 where it affects 
only the unobservable mixing of right-handed quarks. The great difference in magnitude 
between Vcb and sin^ 29f^r is thus a consequence of D^ and L" being highly asymmetric and 
the peculiarities of SU{5) invariance. 

As shown in [^ and [Q, the matrices in Eq. (1) give a remarkably good fit to all the 
known quark and lepton masses and mixings, in particular, when the small-angle MSW 
solution [0] is relevant for the solar neutrino problem. We now review the lepton results and 
explore the possibility of large-angle z/g — z/^ mixing in this 50(10) unified model. 



The lepton mixing matrix is given by 

Vlepton = Uj^U^, (2) 

where Ul is the unitary transformation of the left-handed charged leptons required to di- 
agonahze L°, and 11^ is the complex orthogonal transformation of the left-handed neutrinos 
required to diagonalize the light-neutrino mass matrix, M^ = —N^Mj^^N. In 5*0(10) the 
Dirac neutrino mass matrix A^° and up quark mass matrix f/° are related, and as given 
in Eq. (1) have vanishing first rows and first columns. This is to be regarded as only an 
approximation to the real world, but as we shall see later, it is a very good approximation. 
With this particular texture for A^°, no matter what form Mf> assumes, the light-neutrino 
mass matrix Mj, will also have vanishing first row and column, and we can write 
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where the parameters s and c are complex, in general, with c^ + s'^ = 1. From the form 
of A^° it is easy to see that, formally speaking, \s\ = 0{e). If the unknown matrix M^ is 
parametrized by {M^^)ij = ajjA^\ then, in units of (M^/Aji), one has {Mt,)ij = {Mi^)ji = 0, 
(M^)22 = e^oss, {M^)23 = (M^)32 = eoss -e^a23, and (M^)33 = 033 -2ea23 + 6^022. This gives 

tan 2^2^3 - 2.c/(c^ - s^) = 2e ( /"'T" , ) , (4) 

Va33 - 2ea23 + £^022 - 6^033/ 

and thus ^23 ~ ^ > unless the parameters Uij are fine tuned to have a special relationship to 

each other. 

We see then, that the contributions to the leptonic mixings coming from Ui, are either 

zero or small. The leptonic mixings thus arise from Ul- This is good news, since L° in Eq. 

(1) is known in this model. With the values of the parameters e, a, 6, and 6' given earlier as 

determined by fitting known quantities, one finds that 
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S23/C23 = cr + e = 1.9, 
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Note that the mixing 623 is large and the mixing 6^2 i^ small. As can be seen from Eq. (1), 
the smallness of the angle 6^2 is related to the smallness of the corresponding angle for the 
quarks. The largeness of 62^, on the other hand, is due to the large "lopsided" entry a. There 
is no similar entry for the mixing of the first family. One could imagine introducing one, 
but one would find that doing so would make it hard to fit many known quantities such as 
Vus, Vub, fne/m^, and m^/ms. It would seem, then, that this model must give large z/^ — u^ 
mixing and small z/g — v^ mixing. However, as we shall now show, the same model actually 
can give bimaximal mixing. 

We will suppose now that the up quark matrix element U^i is not exactly zero, but is 
given by tjU^^. If wi„(l GeV) ~ 4 MeV, then 77 ^ 6 x 10~^, which is a thousand times 
smaller than the smallest of the other model parameters, 5 and 5' . It is in this sense that 
the vanishing of the first row and column of U^ is an excellent approximation. In S'0( 10) 
the simplest possibility is that N^^ is also given by rjMu- Thus, we assume for N^ the form 
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This tiny modification of Eq. (1) allows interesting effects for certain forms of Mr. What 
we find is that if the off-diagonal elements in the first row and first column of Mr are 
small or zero, then the small-angle MSW solution results as in p , whereas if these elements 
are important, the large-angle solution of either the "just-so" vacuum or large-angle MSW 
oscillation type can result. Instead of looking at the most general forms for Mr, we will 
illustrate this in two representative cases: 
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Small-angle MSW Solution 

Form (I): 



Mi = -iV°^(M^)-iiVO 





V 



Ml 



A 



R 



(9) 



\ 

-e^A 

t^jA -2tlA - eV^' / 
This liglit neutrino mass matrix is diagonalized by an orthogonal transformation of the type 
given in Eq. (3) for which Eq. (4) applies with ^23 ~ ^- The presence of the new t] 
contribution to the 11 element of A^'' thus does not modify the small- angle MSW solution 
obtained earlier, if M^i is of type (1). The ratio of the mass differences required for this 
solution, 5m\^ = 3.5 x 10"^ eV^, 5m^2 — 7 x 10"^ eV^, can easily be obtained with A ^ — e, 
while the absolute mass scale follows with Mjj ~ 100 GeV and Ar ~ 1.5 x 10^^ GeV. 
Moreover, the mixing results, sin^ 26'^^ = 1.0 and sin^ 26*6^ = 0.008, lie in the desired ranges 
for the atmospheric and small-angle MSW solutions ||T], 0. 

Vacuum Oscillation Solution 

Form (II): 
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A rotation in the 2-3 plane by an angle tan ^23 = ^ will eliminate the 22, 23, and 32 entries 
and induce 12 and 21 entries that are equal to {ri/A)e'^ (neglecting terms higher order in e^). 
Following this with a rotation in the 1-3 plane by an angle O^^ = {ri/A)e brings the matrix 
to the form 
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For \ri/A\ <^ 1, a pseudo-Dirac form for the mass matrix of z/g and u^ obtains with nearly- 
degenerate neutrinos. One finds 
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7 X 10~^, which we shall ignore, and 
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With the use of Eqs. (2) and (5), this gives for the lepton mixing matrix 
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with S23 = sin6'^^, where 6*^^ 

sin^26'e^ — 1.0, which are consistent with the experimental limits. 

It is interesting that the very small value of 5ml2 needed for the vacuum oscillation 
solution to the solar neutrino problem has a natural explanation in this approach. From Eq. 
(12) it is apparent that m^^ = rriu^ oc 1]. In other words, the very tiny parameter rj required 
to fit mu/mt — a quantity pertaining to the first family — actually ends up controlling the 
masses of both the muon- and electron-neutrino which are nearly degenerate. 

Large-angle MSW Solution 

The question arises whether form (II) can also give the large-angle MSW solution to 
the solar neutrino problem 0. The answer is yes, though in this case sin^ 29 ef_i departs 
significantly from maximality. For this solution one needs to have Ar7i^2 — 2 x 10~^ eV^ along 
with Am23 — 3.5 x 10~^ eV^. Although Eq. (11) still applies, the condition for pseudo-Dirac 
neutrinos is no longer satisfied as one finds r]/A ~ 1.8. Three Majorana neutrinos emerge 
for which 
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With Mfj ~ 100 GeV and A^j ~ 1.7 x 10^^ GeV, the three neutrino masses are given 



numerically by 5.9 x 10 ^ eV, 4.7 x 10 ^ eV and 9.3 x 10 '^ eV. Making use of the above 
value for rj/ A, we can write down the equivalent of Eq. (13) for the neutrino mixing matrix 



/ 0.389 -0.878 -1.73e \ 
f/^ = 0.922 0.392 0.961e 
V -0.18e -2e 0.961 / 

The lepton mixing matrix is then found numerically to be 
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( 0.360 -0.908 -0.186 \ 

0.469 0.367 -0.812 
V 0.811 0.222 0.556 / 



(17) 



from which one obtains sin^ 2^^^ = 0.82 and sin^ 2^e/i = 0.43. These values for the mixing 
angles are on the low side of the allowed experimental region for this large-angle MSW type 
of solution. Moreover, the A ~ 3 x 10~^ parameter required is some thirty times smaller than 
that found earlier with a pair of pseudo-Dirac neutrinos; thus a considerably larger hierarchy 
is required in the right-handed Majorana neutrino matrix to reproduce the large-angle MSW 
mixing than for the vacuum solution. These features have also appeared in other forms we 
have assumed for Mr. 

In summary, we have shown that, through the introduction of a small correction which 
gives mass to the up quark and hence also modifies the related Dirac neutrino matrix, 
a bimaximal solution to the solar and atmospheric neutrino oscillations can be achieved, 
provided the right-handed Majorana neutrino matrix mixes the first family with the other 
two. The large-angle vacuum solution is somewhat preferred over the large-angle MSW 
solution for the solar neutrino problem, since a smaller hierarchy is required in the Majorana 
matrix and the mixing angles are more nearly maximal as suggested by present experimental 
data for those type of solutions. But as presently understood, the model does not suggest a 
preference for the large-angle solutions over the small-angle MSW solution. 
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